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Abstract 

EDDE is a Monte Carlo event generator for different Exclusive and Semi-Inclusive 
Double Diffractive processes. The program is based on the extended Regge-eikonal 
approach for "soft" processes. Standard Model and its extensions are used for "hard" 
fusion processes. 
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1 Introduction 



In the present paper we give a brife description of the Monte-Carlo event generator EDDE, 
which includes Exclusive Double Diffractive Events (EDDE, p+p — * p + M +p) and Semi- 
inclusive Double Diffractive Events (SI DDE, p + p ^ p + {XMY} + p). Here M is the 
central particle or system of particles (Higgs boson, RSI model particles, graviton, Xc,b, 
glueball, QQ, gg, 77, QQg, ggg), X and Y denote "soft" gluon radiation in the central 
region, and "-I-" means Large Rapidity Gap. 

EDDE gives us unique experimental possibilities for particle searches and investiga- 
tions of diffraction itself. This is due to several advantages of the process: a) clear 
signature of the process; b) possibility to use "missing mass method" that improve the 
mass resolution; c) background is strongly suppressed; d) spin-parity analysis of the cen- 
tral system can be done; e) interesting measurements concerning the interplay between 
"soft" and "hard" scales are possible [1]. All these properties can be realized in common 
CMS/TOTEM detector measurements at LHC [2]. 

SI DDE is important as a source of main backgrounds for the exclusive processes. 

Large number of event generators are devoted to partonic processes of the Standard 
Model and its extensions, i.e. work at small distances. It is wellknown, that perturbation 
theory has some problems in the description of processes at large distances. That is why 
diffractive processes are usually considered as special cases which description is based on 
different phenomenological approaches. The most popular approach is the Regge-eikonal 
model. 



2 Physics of EDDE 




The exclusive double diffractive process is related to the dominant amplitude of the 
exclusive two-gluon production. Driving mechanism of this processes is the Pomeron. For 
calculation of cross-sections we use the method developed in Refs. [3], [4]. It is based on 
the extension of the Regge-eikonal approach, and succesfuUy used for the description of 
the HERA [5], [6] and p + pip) p + p{p) [7] data. 
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In the framework of this approach, amphtudes for EDDE and SIDDE can be sketched 
as shown in Fig. 1. 

In the model we consider the kinematical region which is typical for diffraction: 



0.001 GeV"^ < |ti2| < 5 GeV^ , 



S 



0.3 



(2) 



where 2 = transfer momenta squared and ^i_2 — ^Pl 1,2/ y/s, ^Pl 1,2 are 

longitudinal momentum transfers of protons. 

The off-shell gluon-proton amplitudes Ti^2 are obtained in the extended unitary ap- 
proach [8] . After the tensor contraction of these amplitudes with the gluon-gluon fusion 
vertex, the full "bare" amplitude Tm depicted in Fig. la inside the dashed oval looks hke 



T, 



M 



TT 
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ap(0) 
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Here 



a^(0)ln 



M 



(3) 

(4) 
(5) 



with the parameters of the "hard" Pomeron trajectory, that appears to be the most 
relevant in our case, presented in the Table 1. The last factor in the r.h.s. of (3) is 



ap{ti)+ap{t2) 



h{v,t) 
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n=2 



exp 



:-i 



nl ■ n \ 
'bi(n-l) 



M bi{v) 
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i7r(ap(0) - 1) 
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,ap{0)-l 



n 



V 

bi 



a'p{0) In^; -|- 60 



(7) 

(8) 
(9) 



Here = —q^ ~ q^, /i = M/2, and sq is a scale parameter of the model which is also 
used in the global fitting of the data on pp (pp) scattering for on-shell amplitudes [7]. 
The fit gives sq 1 GeV^. In (6) we take into account rescattering corrections for "soft" 
amplitudes Ti_2, which play role in the investigations of diffractive pattern [1], [9] for 
t > 1 GeV and give very small contribution to the total cross-section. Below we perform 
also calculations of the so called "soft survival probability", which corresponds to "soft" 
rescattering corrections in the initial and final states. 
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If we take into account the emission of virtual "soft" gluons, while prohibiting the real 
ones, that could fill rapidity gaps, it results in a Sudakov-like suppression [10]: 



Fs{l\i^') = exp 



277 



1-A 



j zPgg{z)dz + jY,Pq9{z)dz 



z{l - z) 



A = 



Pt 



Pt + iJi 



(10) 

(11) 
(12) 



The "hard" part of the EDDE amplitude, Fgg^M, is the usual gluon-gluon fusion 
amplitude calculated perturbatively in the SM or in its extensions. 



Table 1: Phenomenological parameters of the "hard" Pomeron trajectory obtained from 
the fitting of the HERA and Tevatron data (see [3], [6], [9]), and data on pp (pp) scatter- 
ing [7] .The value of the Cgp is corrected in accordance with the latest data from CDF [11]. 



ap{0) 


a'p{0), GeV-2 


rl, GeV-2 


rl, GeV-2 


Cgp 


1.203 


0.094 


2.477 


2.54 


3.2±0.5 



The data on total cross-sections demand unambiguously the Pomeron with larger- 
than-one intercept, thereof the need in unitarization, i.e. taking into account "soft" 
rescattering in the initial and final states. The full amplitude in the Fig. la with unitary 
corrections, T]^'*, is given by the following analytical expressions: 



rpunitar 



(pi,P2, Ai, A2) 



1 f d'^q 
US' J ?2^ 



d^qr d^^rp 



V{s,qT) 



X 



16 ss' J (27r)2 (27r)2 
Tm(pi - qT,P2 + Qt, Ait, ^2t) V{s', q^) , 



V{s,qT) = As{27rYS\qT)+AsJ d^hi 



(13) 
(14) 



where Ai^ = Ai — qt — q'x, A2T = A2 + Qt + q't^ ^^'^ th^ eikonal function Spp^pp can be 
found in Ref . [7] . Left and right parts of the diagram in Fig. la denoted by V represent 
these "soft" re-scattering effects. As was shown in [12], these "outer" unitary corrections 
strongly reduce the value of the corresponding cross-section and change an azimuthal 

angle dependence. 

To calculate differential and total cross-sections for exclusive processes we can use the 
formula 



da 
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(15) 
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M 



Jd'Aid'A2\TM\' 



(17) 



is the "hard" exclusive singlet gluon-gluon fusion cross-section 
"soft" survival probability. For fixed energy is approximately 



where da-^""^^ /d^gg^M 
and 5*^ is the so called 
constant. 

We can extend our approach to the case of additional "soft" radiation in the central 
rapidity region. This process is depicted in Fig. lb. First of all wc have to calculate 
unintegrated gluon distribution fg/g{x,kf, ^'^) inside a gluon. For this task we use the 
method similar to the one presented in [13] (convenient for further Monte-Carlo simula- 
tion) and [14]. Generation is performed step by step by the initial gluon splitting. In 
this case Sudakov-hke form factor Fg in the above formulaes for cross-sections has to be 
replaced by the new one [9] 



i-2 



^x, k^, 



[18) 



and partonic exclusive gg ^ M cross-section has to be replaced by the inclusive one. Since 
there are some theoretical uncertainties concerning the choice of the so called " factoriza- 
tion" scale II, this question has to be considered explicitly by methods of the renormaliza- 
tion theory. In the next version of the generator we will present more exact calculations 
of fg/g{x,k^, fi'^) based on the Bethe-Salpeter equation, taking into account also color 
correlations, which are included in this version in a more simple way. 



3 Main subroutines and common blocks 

Basic functions and subroutines: 

• Subroutine EDDETTPHI(NX,MX,GT1 ,GT2,GFI0) generate h, t2, (p - distribution for 
different masses MX of the central system and different quantum numbers of the 
system: 

NX= 1 — > 0++ (default) 
NX= 2 — ^ = 0, ±2 
NX= 3 — ^ 0"+ 
NX= 4 — ^ glueball 
NX= 5 — > J^ = ±2 

• Function EDDEX(MX) generate ^ distribution for a final proton. 

• Functions of the type DCSxx(N,M,ETA) or DCSQQ(MQ,N,M,ETA) compute exclusive 
differential partonic cross section 

dr]* ' 
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where 7]* is the pseudorapidity of the final parton in the initial gg central mass 
frame. MQ is the mass of a quark. Functions CSQQ(MQ,M) or CSxx(M) are integrated 
exclusive gluon-gluon fusion cross-sections. xx=GG,2GAM,3G,QQG denotes the final 
system {gg,jj,ggg,QQg)- N is an auxiliary integer parameter. If N= then we obtain 
exact cross-section, else functions return an upper estimations for cross-sections. 

• Functions of the type DCSxxSI(N,M, ETA) or DCSQQSI(Mq,N,M,ETA) are inclusive dif- 
ferential partonic cross-section (without the rule = 0). Functions CSQQSI(MQ,M) 
or CSxxSI (M) are integrated inclusive gluon-gluon fusion cross-sections. 

• Subroutines of the type GENERyyxx(GMGG,GETAJ) , GENERyyQQCMQ.GMGG.GETAJ), 
GENEREX3G(GMGG,GETAJ,GX3,GPT3,GFI3) {ggg distributions), 

GENEREXQQGCGMGG , GETAJ , GX3 , GFIS , GTHETAS) {QQg distributions) generate distri- 
butions in 7^* (GETAJ) and mass of the central system (GMGG). Also GX3, GPT3, GFI3, 
GFIS, GTHETAS are variables for the 3rd jet in the 3-jet kinematics. yy=EX, SI means 
Exclusive and Semi-inclusive events. xx=GG, 2GAM {gg and 77 final systems). 

• Subroutine SICASCAD2(MC,N1,PG1,P1,N2,PG2,P2,MX,NFAIL) generates two corre- 
lated "soft" gluonic systems X and Y in the process p + p ^ p + {XMY} +p. MC is 

the mass of the final system (Higgs boson, Radion,graviton,glueball,Xc,6,5'5',QQ, 77), 
MX is the mass of the system {XMY}, Nl, N2 are numbers of gluons in X and 
Y. PG1(5,500), PG2(5,500) are four momenta of gluons in X and Y, and PI (5), 
P2(5) are momenta of "hard" gluons in the process gg — > M. 

• Function SOFTSURV(N,MGG) returns the "soft survival probability" factor S"^. MGG is 
the central mass. 

• Function M2DLUMDM2(N1,M2, MGG) returns "luminosity" L^^^^- ^^^^^ integrated in 
the rapidity of the central system. Nl= for EDDE, Nl= 1 for SIDDE. N2= for 
resonances, N2= 1 for jets and 77. 

• Subroutine EDDERS1C(NRS,RSXI,RSGAM,RSMH,RSMR,RSM0BS,RSWD, BR) returns mas- 
ses, widths and branching fractions for resonances of the RSI model [15]. Parameter 
NRS switches between two mass states H* (Higgs boson, NRS= 1) and R* (Radion, 
NRS= 2). RSXI, RSGAM, RSMH, RSMR are input parameters for the RSI model (mixing 
parameter ^, scale parameter 7, "bare" masses m^, rrir before mixing) [15]. RSMOBS, 
RSWD, BR are the mass, width and branching fraction into bb state for H* or R*. 

• Function EDDECS(IP) retuns total cross-section of the corresponding process: 

IP= 400: EDDE, resonance production (Standard Model Higgs boson, H* and R* of RSI 
model, Xc,b, glueball, graviton). Partonic gg ^ H cross-section could be found 
in [16]. Other EDDE cross-sections were calculated in [3], [17]. 

IP= 401: EDDE, QQ production. 

IP= 402: EDDE, gg production. 

IP= 403: EDDE, 77 production. 
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IP= 


404 


EDDE, QQg production. 


IP= 


405 


EDDE, ggg production. 


IP= 


406 


SIDDE, resonance production. 


TP— 


407 


STDDE no nrndiiptinn 


IP= 


408 


SIDDE, gg production. 


IP= 


409 


SIDDE, 77 production. 



• Subroutine EDDEEVE returns an event in PYTHIA [18] format according to the value 
of the process identification variable (see data card below) . 

• Subroutine EDDEINI makes initialization of the generator and provide the interface 
with PYTHIA. 

• Subroutine EDDEPUTDAT fills all the internal common blocks. It is called from 
EDDEINI. 

All the cross-sections are calculated in mb. 
EDDE common blocks: 

• Fundamental constants: 

COMMON/EDDEFUND/ MNI , REI , PI , CSMB , LAMQCD , TF , CP , BFO , BFl , NF , NC , NLOSW 

DATA PI/3.141592654D0/, CSMB/0 . 38D+00/ 
DATA REI/ ( 1 . DO , . DO) / , MNI/ (0 . DO , 1 . DO) / 
DATA NC/3/ , NF/5/ ,NL0SW/1/ , LAMQCD/0 . 5D-01/ 
DATA TF/0.5D0/ 

TF,CF,BF0,BF1,NF,NC are SU(3) constants 

NLOSW is a switch parameter. For zero value as is calculated in the Leading Order, 
for unit value we have NLO calculations. 

• Parameters for the Regge-eikonal model and Regge kinematics: 

COMMON/EDDESOFT/ CP (3) , DP (3) , HP (3) , KG (3) , AP (3) , 
& TIMIN , TIMAX , T2MIN , T2MAX , FKK , CGP , NAPR , NFI 

Parameters of trajectories: couphngs CP, intercepts DP= q;/p(0) — 1, RP= r^^, 
RG= rg, slopes AP= a'jp{0) 



DATA CP/0 . 5300D+02 , . 9700D+01 , . 1670D+01/ 
DATA DP/0 . 5800D-01 , . 1670D+00 , . 2030D+00/ 
DATA RP/0 . 6300D+01 , . 3100D+01 , . 2480D+01/ 
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DATA RG/0 . 6300D+01 , . 3100D+01 , . 2540D+01/ 
DATA AP/0 . 5600D+00 , . 2730D+00 , . 9400D-01/ 

ti^2 limits: 

DATA TlMIN/0. lD-02/, TlMAX/0 . 7D+01/ 
DATA T2MIN/0.1D-02/, T2MAX/0 . 7D+01/ 

Number of terms in the expansion: 

DATA NAPR/9/ 
Main constant of gp — > gp amplitude: 

DATA CGP/0.316D+01/ ! (3.2+-0.5) 

NFI=IPAR 5 from the data card, it is the code for azymuthal distribution. 

• Kinematical limits for ^1^2: 

COMMON/EDDETOT/ XI IMIN , XI2MIN , XI IMAX , XI2MAX 

DATA XIlMIN/0 . lD-04/ , XI2MIN/0 . lD-04/ 
DATA XI lMAX/0 . IDO/ , XI2MAX/0 . IDO/ 

• Parameters for "hard" gluon-gluon fusion cross-sections: 

COMMON/EDDEHARD/ MGGCUT , ET JCUT , MXMAX , 
& ETA JMAX , PLUM , PSURV , PSUD , ETASIMAX , SQS , 
& PSIDD1,PSIDD2 

MGGCUT= 2*ETJCUT, ETJCUT=RPAR 2 is the cut on rriT = + m^^rf for final 

"hard" parton. For gluons and photons rripart = 0, i.e. ETJCUT is the cut on the 
transverse momenta. 

MXMAX=500 . DO GeV is the upper bound for the mass of the central system in the 
generator. 

ETAJMAX= max \ri*^rt\ the upper bound for the pseudorapidity of final "hard" par- 
tons in the central mass frame of initial gluons. 

ETASIMAX=RPAR 9 is the pseudorapidity interval for "soft" radiation. 
SQS= ^/s =RPAR 1. 

Other parameters are used for different auxiliary distributions. 

• Auxiliary parameters for 3-jet kinematics: 
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C0MM0N/EDDE3 JP/ DER3 J , XMAX3 J , PAR3G (5) 



• Parameters for RSI model [15]: 

COMMON/EDDERSl/ RSXIO , RSGAMO , RSMHO , RSMRO , NRSO 

RSXIO=RPAR 5= ^ is the mixing parameter, 

RSGAMO=RPAR 6= 7 = 246 GcV/A<^ is the scale parameter, 

RSMHO=RPAR 7= ruh and RSMRO=RPAR 8= rrir are "bare" masses of the model, 

NRSO=IPAR 6 is the switch parameter between SM Higgs boson (IPAR 6= 0), H* 

(IPAR 6= 1), R* (IPAR 6= 2). 

• Additional global parameters: 

COMMON/EDDEOTHER/ KCP , IPROC , AMO , AMP , S , MQ 

KCP=IPAR 4 is the code of the central particle in the process number 400 and 406, 
IPROC=IPAR 1 is the number of process (see above the parameter IP of the function 
EDDECS), 

AMO=RPAR 3 is the mass of the central particle, 
AMP is the proton mass, S=SQS*SQS, 

MQ=RPAR 10 is the mass of heavy quark (b-quark by default). 



• Data tables for luminocities and distribution functions: 

COMMON/EDDETABl/ LUMl (480) ,FLUM1 (30 , 16) , 
& DX1,DY1,X01,Y01 

C0MM0N/EDDETAB2/ RDI3G(630) ,FRDI3G(30,21) , 
& DX2,DY2,X02,Y02 

/EDDETAB3/ RI3GA(480) ,RI3GB(480) , 
& FRI3GA(30,16) ,FRI3GB(30, 16) , DX3 , DY3 , X03 , Y03 

C0MM0N/EDDETAB4/ FT4(10201) ,FMT4(101 , 101) , 
& DX4,DY4,X04,Y04 

C0MM0N/EDDETAB5/ DFT5(10201) ,DFMT5 (101 , 101) , 
& DX5,DY5,X05,Y05 
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4 Data card 

User can manage the process of generation with the help of the data card. Parameters 
that could be changed and their default values: 

• IPAR 1=403: identification number of the process (IP in the function EDDECS). 

• IPAR 2=10000: number of generated events. 

• IPAR 4=25: code of the central particle for processes number 400 and 406 (according 
to PYTHIA codes for particles, for example, 25 for Higgs boson and RSI particles, 
50 for glueball etc). 

• IPAR 5=1: parameter to switch between different azymuthal distributions (NX in the 
subroutine EDDETTPHI). 

• IPAR 6=0: number of particle for the process 400 or 406 and IPAR 4=25. IPAR 6=0: 

Standard Model Higgs boson; IPAR 6=1: H* ("Higgs boson") mass state of the RSI 
model; IPAR 6=2: R* (Radion) mass state of the RSI model. 

• RPAR 1=14000 . OE+00: (in GeV) center mass energy V^- 

• RPAR 2=25 . OE+00: (in GeV) transverse mass cut my = \JE^ + mlart for the "hard" 
final parton of the central system. For gluons and photons this parameter is equal 
to Et cut. 

• RPAR 3=120. OE+00: (in GeV) mass of the central particle for processes number 400 
and 406. 

• RPAR 5=0 . 16E+00: mixing parameter ^ of the RSI model, typical values —0.2 —>■ 0.2. 

• RPAR 6=0.246E+00: scale parameter 7 = 246 GeV/A<^, where A<^ has typical values 
1000 ^ 5000 GeV. 

• RPAR 7=150. OE+00: (in GeV) "bare" mass of Higgs boson before mixing rrih. Typ- 
ical values 120 180 GeV. 

• RPAR 8=1 10. OE+00: (in GeV) "bare" mass of Radion before mixing m^. Typical 
values ~ 100 ^ 300 GeV. 

• RPAR 9=5. OE+00: pseudorapidity interval for "soft" radiation X and Y in the 
SIDDE, -RPAR 9/2 < 7/ <RPAR 9/2. 

• RPAR 10=4.8E+00: (in GeV) mass of the final "hard" quark in the central system. 
Some PYTHIA definitions: 
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PYTHIA 
C 

MRPY 1=77123456 
MSEL 

C . . . some PYTHIA 
MSTP 61 =1 
MSTP 71 =1 
MSTP 81 =1 
MSTP 111=1 
MSTP 122=0 
C . . .Higgs decay 
MDME 210,1 =0 
MDME 211,1 =0 
MDME 212,1 =0 
MDME 213,1 =0 
MDME 214,1 =1 
MDME 215,1 =0 
MDME 216,1 =-1 
MDME 217,1 =-1 
MDME 218,1 =0 
MDME 219,1 =0 
MDME 220,1 =0 
MDME 221,1 =-1 
MDME 222,1 =0 
MDME 223,1 =0 
MDME 224,1 =0 
MDME 225,1 =0 
MDME 226,1 =0 



! State of random number generator 
!full user control 
definitions . . . 

Initial-state QCD and QED radiation 
Final-state QCD and QED radiation 
multiple interaction 
fragmentation and decay 
switch off X section print out 
definition. . . 
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gg 
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-> gamma gamma 
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-> 


gamma ZO 
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-> 


ZO ZO 


hO 


-> 


w+w- 



5 Results from EDDEv2.1 

Here some samples from the work of the generator with PYTHIA [18] are presented ( 
Figs. 2a)-d).) 

The updated version of the generator EDDEv2.1 will be available on the web-page 
http : // sirius.ihep.su/cms/higgsdif f / dif f .html 
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d) 

Figure 2: Some distributions from EDDE. (90 signal and 660 backround events. MTcut = 
25 GeV). Solid curve represent signal+background, dashed one is the background, and 
dash-dotted is the signal, a) azimuthal angle distribution for Higgs boson production; b) 
rapidity and c) pseudorapidity distributions; d) integrated t-distributions. 
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